The K (7T7r)/ =2 Decay Amplitude from Lattice QCD 
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We report on the first realistic ab initio calculation of a hadronic weak decay, that of the amplitude 
Ai for a kaon to decay into two 7r-mesons with isospin 2. We find ReA2 = (1.436 ± 0.063 s ^ a ^ ± 
0.258 S y S t) 10 -8 GeV in good agreement with the experimental result and for the hitherto unknown 

imaginary part we find Imyi2 = — (6.83±0.51 s ^ a ^±1.30 S y S ^) 1CP 13 GeV. Moreover combining our 
result for Im^ with experimental values of Re^, ReAo and e'/e, we obtain the following value for 
the unknown ratio Im Ao/Ke Ao within the Standard Model: ImAo/Re Ao = — 1.63(19) s tat (20) sys t x 
1(T 4 . One consequence of these results is that the contribution from Im^ to the direct CP 
violation parameter t (the so-called Electroweak Penguin, EWP, contribution) is Re(e'/e)EWP = 
— (6.52 ± 0.49 s t a t ± l-24 S y S t) x 10 -4 . We explain why this calculation of A2 represents a major 
milestone for lattice QCD and discuss the exciting prospects for a full quantitative understanding 
of CP-violation in kaon decays. 
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Introduction 



CP-violation is a necessary ingredient for the genera- 
tion of the matter-antimatter asymmetry in the universe 
and understanding its origin, both within and beyond the 
Standard Model, is one of the primary goals of particle 
physics research. It was first discovered in the decays of 
kaons into two pions, and in this paper we report on a 
calculation of the amplitude A2 for K — > (1111)1=2 decays 
from first principles. (Bose symmetry implies that the 
two-pion eigenstates have isospin or 2 and we denote 
the corresponding complex amplitudes by Aq and A%.) 
For KeA2 we find good agreement with the known ex- 
perimental value (see Eq. (jlOp l and, more importantly, 
we are also able to determine the previously unknown 
quantity ImA2 (Eq. (jll[) l . A2 is obtained by combining 
our lattice results for matrix elements of the four-quark 
operators in the effective Hamiltonian (the matrix ele- 
ments are given in Eqs. ©-(JSl) ) with Wilson coefficients 
and CKM-matrix elements. In addition, within the Stan- 
dard Model we can combine our result for Im A2 with 
the experimental values of Rev4 , Re A 2 and e'/e to de- 
termine the remaining unknown quantity IthAq, so that 
both the complex amplitudes Aq and A 2 are now known. 

This is the first quantitative determination of an am- 
plitude for a realistic hadronic weak decay and extends 
the framework of lattice QCD into the important domain 
of non-leptonic weak decays. It has taken several decades 



for a realistic lattice QCD calculation of K — > (irir)j = 2 
decay amplitudes to become possible because very sig- 
nificant theoretical developments and technical progress 
were required. These are briefly discussed below and 
explained in more detail in |1J, where a full description 
of our calculation can be found. Among the key issues 
is the fact that performing the simulations in Euclidean 
space makes the evaluation of 7T7t rescattering effects non- 
trivial. The presence of two pions interacting strongly in 
a finite box leads to finite-volume effects which must be 
controlled Q • In fact the effects of finite volume can be 
exploited to ensure the equality of the energies of the ini- 
tial kaon and final two-pion states by the imposition of 
carefully devised boundary conditions. Finally, it is only 
relatively recently, with the improvement of algorithms 
and access to teraflops-scale computing resources, that it 
has become possible to perform simulations at physical 
u and d-quark masses. 

The calculation of A2 is also important in that it de- 
termines the 0(5%) contribution of direct CP violation 
to e 0, 0] , a level of precision which has become relevant 
due to the major recent improvements in the evaluation 
of the Bk parameter with an uncertainty of less than 3%, 
see e.g. [i[ (for a recent review see @). Of course, a com- 
plete understanding of CP-violation in K — > irir decays, 
including the evaluation of e'/e and an understanding of 
the AI =1/2 rule, requires the ability to compute both 
Aq and A2 directly. The present calculation is an im- 
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portant milestone on the road to achieving this. For Aq 
however, the two pions have vacuum quantum numbers 
and the correlation functions are dominated by the vac- 
uum intermediate state. We report on our exploratory 
work in developing techniques for the efficient subtrac- 
tion of the vacuum contributions in Q and look forward 
to presenting results of a realistic computation of Aq in 
the future. 



Details of the Simulation 

The simulations were performed using domain wall 
fermions (DWF) with a gauge action which we call 
IDSDR. This is the Iwasaki action modified by a weight- 
ing factor, called the Dislocation Suppressing Determi- 
nant Ratio (DSDR) gHIoj], which allows us to suppress 
configurations with large numbers of modes of the 5- 
dimensional DWF transfer matrix with near unit eigen- 
value while retaining adequate topological change. This 
modification is necessary since we have a relatively large 
lattice spacing a which increases the frequency of dislo- 
cations which break the chiral symmetry. 

We have generated two ensembles of 2 + 1 flavor DWF 
with the IDSDR gauge action at f3 = 1.75 and a lat- 
tice size of 32 3 x 64 x 32, where the final number is the 
length of the fifth dimension. We determine the resid- 
ual mass to be m ros = 0.00184(1) [ll|. (Masses writ- 
ten without units are to be understood as being in lat- 
tice units.) The ensembles are generated with a simu- 
lated strange-quark mass of = 0.045 and light-quark 
masses of m; = 0.001 and mi — 0.0042, with correspond- 
ing unitary pion masses of approximately 170MeV and 
250 MeV respectively. Quark propagators are generated 
for a range of valence masses. The analysis presented in 
this paper is performed using 63 configurations from the 
0.001 ensemble, each separated by 8 molecular dynam- 
ics time units, and quark propagators with = 0.049 
and mi — 0.0001. The (partially-quenched) pion has a 
near-physical mass of approximately 140 MeV. A subse- 
quent detailed analysis with greater statistics and im- 
proved procedures have yielded a slightly lower value for 
the bare physical strange quark mass (0.0464(7)) [Tlj | . 

We obtain the lattice spacing and the two physi- 
cal quark masses m uc i and m s using a combined anal- 
ysis of these IDSDR ensembles and our (3 — 2.25, 
32 3 x 64 x 16 and /3 = 2.13, 24 3 x 64 x 16 domain 
wall fermion configurations with the Iwasaki gauge ac- 
tion [H, EH • This involves a combined fit of pion and 
kaon masses and decay constants and the mass of the 
fi-baryon as functions of the quark masses and lattice 
spacing. We extrapolate to the continuum limit along 
a family of scaling trajectories defined by constant val- 
ues of m Vl mx and m^ fl3j ] - In our fits we take the 
lattice artefacts to be 0(a 2 ) as expected and note that 
the coefficients of the a 2 terms are not equal for the 



two different lattice actions. From the combined chi- 
ral and continuum fits we obtain for the IDSDR ensem- 
bles a -1 = 1.375(9) GeV and physical quark masses of 
mi = 0.00174(3) and m s = 0.0483(7) in lattice units, 
which correspond to m^f (2 GeV) = 3.43(13) MeV and 



(2 GeV)=95. 1(3.0) MeV. (Here m = m + m res .) 



In order to ensure that the energy of the two-pion fi- 
nal state (in the rest-frame of the kaon) is (almost) equal 
to rn_R- , we have carefully chosen both the volume of the 
lattice and the boundary conditions on the quark fields. 
With periodic boundary conditions, the two-pion ground 
state corresponds to each pion being at rest (up to finite- 
volume corrections) so for the physical decay we would 
need to consider an excited state [2|. Instead we intro- 
duce antiperiodic spatial boundary conditions for some 
components of the e?-quark's momentum, so that the cor- 
responding components of the momentum of a n + meson 
are odd- integer multiples of ir/L (L = 32 is the spatial 
extent of the lattice) . There is now no state with both 
pions at rest. This is not sufficient however, since the 
physical decay K + — > ir + n° involves a 7r° which, even 
with antiperiodic boundary conditions on the d-quark, 
has momentum components which are integer multiples 
of 2tt/L. It is therefore not possible to construct the 
7T + 7r° state at rest. These problems are overcome by us- 
ing isospin symmetry and the Wigner-Eckart theorem to 
relate the matrix elements for the decay K + — > 7r + 7r° to 
those for the unphysical process K + — > ir + n + : 
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= 1/21-"- I 2 V l v *A/ a =3/2l 

an approach proposed and first explored in 14 , 15[ . The 
superscripts and subscripts on the operators Q denote 
how the total isospin I and the z-component I z change 
between initial and final state. Neglecting violations of 
isospin symmetry, ([T]) is exact and so we can use the 
two-7r + state to calculate the physical AI — 3/2 decay 
amplitudes. In this way we are able to avoid the need to 
consider an excited state and at the same time we reduce 
the required size of the lattice. 

The kaon and pion masses and the energy of the two- 
pion state in the simulation, together with the corre- 
sponding physical values are presented in Tab. [J 

Evaluation of Ai 

The generic form of the effective Hamiltonian for K — > 
7T7T decays is 



G 



H c ff = —j= ^2 (^ckm)i Ci Qi , 



(2) 



where Gf is the Fermi constant, (Vckm)i are the ap- 
propriate CKM-matrix elements, (specifically we use 
V us = 0.2253, V ud = 0.97429 and r = -V t * s V td /V* s V ud = 
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m K + m w+ 


Etttt ma — E-kk 


Simulated 
Physical 


511.3(3.9) 142.9(1.1) 
493.677(0.016) 139.57018(0.00035) 


492.6(5.5) 18.7(4.8) 

m K+ 



TABLE I: m K + , and E^-n in the simulation and the corresponding physical values. The results are given in MeV. 



0.0014606 - 0.00060408i), Q t are four-quark operators 
and Ci are the Wilson coefficients. The calculation of 
A2 requires the evaluation of the matrix elements of 
three operators, classified by their transformations un- 
der SU(3)l x SU(3)_r chiral symmetry: 

0(274) = )l (u j di) L , Q m , = (?tT) L (u j di) R , 
Q(8 < 8) m ^=(s i d j ) L (u^d i ) R , (3) 

where i, j are color labels which run from 1 to 3. (Qrs,e) 
and Q(8,8)mix are the EWP operators contributing mainly 
to Im^-) The main achievement being reported here 
is the successful determination of the matrix elements 
i=2{^^\Qi\K). This starts with the evaluation of the 
correlation function 

C K ^(t K ,t Q ,t^) = {{)\J^{t^)Q l {t Q )J K {t K )\Q) 

= e -m K {t Q -t K ) e -E„„(t„„-t Q )^Q j j^q) 1^ x 

(7nr\Q i (0)\K)(K\4(0)\0)+--- (4) 

where J K and J T7r are interpolating operators for the 
kaon and two-pion states, which are summed over space 
and hence have zero momentum. The energy of the two- 
pion state, E^, is a little larger than I^Jvrv^ + n{ir/L) 2 
because of finite- volume effects (in the isospin 2 state the 
two-pion potential is repulsive). Here n is the number of 
spatial directions in which anti-periodic boundary con- 
ditions have been imposed on the d-quark. The ellipses 
represent the contributions of heavier states, which are 
suppressed if tQ — tx and t^ n — tQ are sufficiently large. 
The sources for the kaon and two-pions are placed at 
fixed times, tx and t n7r (in lattice units), and we vary 
the position of the operator tQ. 

The required (7r7r|Qj|if) matrix element is one of the 
factors in Eq. @ and we need to remove the remain- 
ing factors. This is achieved by evaluating two-point 
correlation functions Cx(t) = (0 | Jjr(t) jt-(O) | 0) and 
C n7r (t) — (0 I Jlni®) I 0), and calculating the ratio 

R(tC)) = Cxmrjtx, tQ, t^n) ^ 

^ (™\Qi\K) 

(0| J^(0)\wir)( K I J K (0) I 0) ' 

where the factors in the denominator of Eq. ^ are deter- 
mined by fitting the correlation functions Ck and C^. 
R(tQ) is independent of tQ if all the time intervals are 
sufficiently large. For illustration of the plateaus we 



present in Fig. Q]the tQ behavior for the 3 operators for 
t^jr — tx = 24 . (We also have results for t^ = 20, 28 
and 32.) 

Having obtained the matrix elements of the bare lattice 
operators ( / Kir\Q\ a " tt \K) , in order to obtain A 2 we must 
rcnormalize the operators and apply finite- volume correc- 
tions. The latter are given by the Lellouch-Liischer factor 
in terms of the s-wave Trrr-phase shift |2| (the phase-shift 
can be obtained from E n7r 16]). In order to combine 
our results with the Wilson coefficients calculated in the 
MS-NDR scheme (l7l - [l9j |. we perform the renormaliza- 
tion in 3 steps. We start by obtaining the renormaliza- 
tion constants in four RI-SMOM schemes using the pro- 
cedures described in |5| . Because the lattice is coarse the 
renormalization scale is chosen to be low, 1.145 GeV, to 
avoid lattice artefacts. We determine the universal, non- 
perturbative continuum step scaling function required to 
evolve the operators to 3 GeV using our Iwasaki lat- 
tices [20|, [2l| . Finally at 3 GeV we convert the results 
to the MS-NDR scheme using one-loop perturbation the- 
ory. 

Our final results for the matrix elements in the MS- 
NDR scheme at a renormalization scale of 3 GeV are: 

M (27 ,i) = (3.20 ± 0.13 sta t ± 0.58 syst ) 10~ 2 GeV 3 , (7) 

M m = (5.85 ± 0.89 sta t ± 1-llgyst) 10" 1 GeV 3 , (8) 

M (8 , 8)mix = (2.75 ±0.12 stat ±0.52 syst ) GeV 3 , (9) 

where for each operator Qi, Mi = (tt + tt + \ Qi\K + ) . 
In terms of these matrix elements, A 2 e % 2 = 
v / 3/2(G F /v / 2) Ei (VcKM)i Q M h where the Wilson co- 
efficients correspond to operators for the physical K + — > 
7r + 7T° decays with the normalization (sg?)l[(uu)l — 
(dd)]L + (su)i(ud)i for the (27, 1) operator and similarly 
for the EWP operators. 

Combining the results in Eqs. (jTj) - (J9J) with the Wilson 
coefficients, CKM matrix elements and Gf we find: 

Re A 2 = (1.436±0.062 sta t±0.258 syst ) 10" 8 GeV (10) 
lmA 2 =-(6.83±0.51 stat ±1.30 syst )10- 13 GeV. (11) 

The result for Re A 2 agrees well with the experimental 
value of 1.479(4) x 10~ 8 GeV obtained from K + decays 
and 1.573(57) x 10~ 8 GeV obtained from K s decays (the 
small difference arises from the unequal u and d quark 
masses and from electromagnctism, two small effects not 
included in our calculation). Im^ is unknown so that 
the result in Eq. (fTTj) provides its first direct determi- 
nation. For the phase of A 2 we find lm. A 2 /~ReA 2 = 
-4.76(37) stat (81) syst 10- 5 . 
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FIG. 1: R(tQ) for the 3 operators which contribute to K — !> (7T7r)j-=2 decay amplitudes: (a) (sc()_l (u<1)l, (b) (sd)i (ud)n and 
(c) (s'cP)/. (uid l )n, where (sg()l (ucI)l,r = (s l 7 M (l — 7 5 )d*) (u-'7 M (l =F7 )d J ). i,j are color labels and tx and t^x are and 24. 



The various sources of systematic error are analysed 
in detail in [1] and our conclusions are summarised in 
Tab. UIl The dominant source of uncertainty is due to 
lattice artefacts, and since we have a relatively coarse 
lattice and the matrix elements are proportional to a~ 3 , 
these errors are substantial. The estimate of 15% is ob- 
tained in two ways: from the variation in the value of 
a obtained using itiq, f^, Jk and r to set the scale 
and from the a 2 term in global chiral-continuum fits of 
the Bk parameter of neutral kaon mixing (fits are per- 
formed using both IDSDR and Iwasaki lattices). The 
finite- volume uncertainties are estimated from the differ- 
ences of infinite- and finite-volume one-loop chiral per- 
turbation theory. The uncertainties in the Wilson coeffi- 
cients are conservatively taken as the difference between 
the leading and next-to-leading order terms as defined 
in 22j. We estimate the truncation errors in the pertur- 



bative factors converting the operators to the MS-NDR 
scheme from the variation of the results obtained using 
different RI-SMOM intermediate schemes. We note also, 
that in contrast to AI =1/2 decays, all the quarks par- 
ticipating directly in AI — 3/2 decays are valence quarks 
and in such cases the effect of using partially quenched 
or partially twisted boundary conditions are small 23| . 



For more details and for a discussion of the remaining 
uncertainties, due to the small difference from physical 
kinematics, and in the evaluation of the Lellouch-Liischer 
factor and the step-scaling functions, we refer the reader 
to d. 

Our result for Im A 2 can be combined with the experi- 
mental results for Re A 2 , Re A Q = 3.3201(18) x KT 7 GeV 
and e'/e to obtain the unknown ratio: 



^ = -l-63(19) stat (20) 
Re Aq 



syst 



10" 



(12) 



This ratio allows us to determine in full QCD the ef- 
fect of direct CP violation in Kl — > tttt on e, customarily 
denoted by K e 3], (K e ) abs = 0.923±0.006. where the sub- 
script "abs" denotes that at present only the absorptive 
long-distance contribution (Im T12) is included Q (the 





ReA 2 


lmA 2 


lattice artefacts 


15% 


15% 


finite- volume corrections 


6.2% 


6.8% 


partial quenching 


3.5% 


1.7% 


renormalization 


1.7% 


4.7% 


unphysical kinematics 


3.0% 


0.22% 


derivative of the phase shift 


0.32% 


0.32% 


Wilson coefficients 


7.1% 


8.1% 


Total 


18% 


19% 



TABLE II: Systematic error budget for Reyl2 and Im42 



error is now dominated by the experimental uncertainty 
in e'/e). The analogous contribution from the disper- 
sive part (Im M12) [4J is yet to be determined in lattice 
QCD, but we describe progress towards being able to do 
this in 

Using our value of Im A2 in Eq. (fTTj) and taking the ex- 
perimental value given above for Re A2 from K + decays 
we obtain the EWP contribution to e'/e, Re(e'/e)EWP = 
-(6.52±0.49 stat ±1.24 syst )xl0- 4 . 



Conclusions and Outlook 

The ab initio calculation of the complex K — > (7T7r)/ = 2 
decay amplitude A 2 described above builds upon sub- 
stantial theoretical advances, achieved over many years 
as outlined in the introduction. It is encouraging that the 
value we find for Re A 2 is in good agreement with experi- 
ment and we are also able to determine Im A 2 for the first 
time. It will be important to repeat this calculation us- 
ing a second lattice spacing so that a continuum extrap- 
olation can be performed thus eliminating the dominant 
contribution to the error, reducing the total uncertainty 
to about 5%. We expect that the dominant remaining 
errors in A 2 will then come from the omission of electro- 
magnetic and other isospin breaking mixing between the 
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large amplitude A and A 2 . 

Much more challenging but of even greater interest is 
the application of these methods to the evaluation of A 
allowing for a calculation of e'/e and an understanding 
of the Al = 1/2 rule. Although the framework pre- 
sented here will also support the calculation of A , se- 
rious obstacles must be overcome. Much larger Monte 
Carlo samples will be required to remove the large fluc- 
tuations remaining after the contribution of the vacuum 
state has been removed. The anti-periodic boundary con- 
ditions for the e?-quark field used in this paper cannot be 
applied to the / — irir state. Instead more sophisticated 
boundary conditions, mixing quarks and anti-quarks and 
an isospin rotation, ( G-parity boundary conditions) [3] , 
must be used for both the valence and the sea quarks. 
Exploratory studies suggest that obtaining adequate 
statistics will be practical with the next generation of ma- 
chines which will become available to our collaboration 
within the next few months. 
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